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Abstract:  Minimum  achievable  temperature  of  ~1 10  K  is  measured  in  a  5% 
doped  Yb:YLF  crystal  at  ~k  =  1020  nm,  corresponding  to  E4-E5  resonance 
of  Stark  manifold.  This  measurement  is  in  excellent  agreement  with  the 
laser  cooling  model  and  was  made  possible  by  employing  a  novel  and 
sensitive  implementation  of  differential  luminescence  thermometry  using 
balanced  photo-detectors. 
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1.  Introduction 


Optical  refrigeration  or  laser  cooling  of  solids  is  based  on  anti-Stokes  fluorescence  [1].  When 
the  laser  excitation  wavelength  is  above  mean  emission  wavelength  (Xf)  of  a  given  transition, 
the  subsequent  fluorescence  requires  phonon  absorption  in  order  to  establish  quasi 
equilibrium.  The  efficient  escape  of  the  upconverted  fluorescence  would  then  carry  heat  out  of 
the  material  resulting  in  net  cooling  [2-4]. 

Laser  cooling  of  a  solid  was  first  demonstrated  in  a  Yb-doped  glass  by  Epstein  et.  al.  in 
1995  [5].  Soon  after  this  report,  many  innovative  applications  such  as  an  all-solid-state  optical 
cryocooler  [6,7]  and  radiationally-balanced  laser  [8]  were  proposed.  To  date,  laser  cooling  has 
been  observed  in  a  wide  variety  of  glass  and  crystalline  hosts,  doped  with  trivalent  rare-earth 
ions  of  Yb,  Tm  and  Er  (for  a  recent  broad  literature  review  see  for  e.g  [3,9,10]).  Excitation 
near  the  Stark-manifold  E4-E5  resonance  of  a  Yb31 -doped  LiYF4  (Yb:YLF)  crystal  has 
recently  led  to  bulk  cooling  from  room  temperature  to  155  K  [11].  Subsequently,  the  same 
system  was  used  to  cool  a  semiconductor  (GaAs)  payload  to  165  K  [11,12],  The  above  results 
were  limited  by  low  pump  power  and  more  importantly  by  the  shortest  wavelength  of  the 
tunable  pump  laser  (Yb:YAG,  1023-1040  nm).  Our  laser  cooling  model  predicted  that  the 
minimum  achievable  temperature  (MAT)  of  ~115  K  should  be  possible  if  the  excitation  is 
tuned  exactly  to  the  E4-E5  Stark  resonance  of  Yb3+  at  -1020  nm  [11].  In  this  paper  we  report 
on  the  experimental  validation  of  the  model  by  conducting  spectroscopic  measurements  of  the 
cooling  efficiency  of  the  Yb:YLF  sample.  The  verification  of  local  cooling  to  110  ±  5  K,  as 
reported  in  Section  4  of  this  paper,  represents  a  major  milestone  in  optical  refrigeration  as  it 
demonstrates  cooling  below  123K  (-150  C),  designated  by  NIST  as  the  onset  of  cryogenics. 
These  measurements  were  made  possible  utilizing  a  novel  and  sensitive  pump-probe 
technique  of  two-band  differential  spectral  metrology  (2B-DSM),  which  probes  local 
temperature  change  directly.  This  is  in  contrast  to  a  more  popular  photo-thermal  deflection 
method  where  various  competing  contributions  complicate  signal  interpretation. 

2.  Modeling  of  the  cooling  efficiency 

The  cooling  efficiency,  defined  as  the  ratio  of  the  cooling  power  ( Pcooi)  to  the  absorbed  laser 
power  {Pabs),  is  given  by  [2]: 


nMJ) 


1 

1  +aJa(X,T) 


Af(T) 


-1, 


(1) 


where  i)ext  is  the  external  quantum  efficiency  (EQE)  defined  as  the  fraction  of  excited  ions  that 
lead  to  a  fluorescence  photon  exiting  the  host  material.  Fligh  rjext  (>  99%)  is  known  for  rare- 
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earth  ions  in  low  phonon  energy  hosts  (such  as  fluorides),  making  them  ideal  candidates  for 
laser  cooling.  Furthermore,  it  is  generally  expected  (and  experimentally  verified)  that  t/ext  is 
rather  temperature  independent  [13].  The  cooling  transition  is  characterized  by  an  absorption 
coefficient  a(k,T)  and  mean  fluorescence  wavelength  kf  (T)  (taking  into  account  effect  of 
luminescence  trapping  and  re-absorption),  while  ab  describes  the  background  parasitic 
absorption.  In  the  sign  convention  adopted  here,  negative  values  of  the  cooling  efficiency 
correspond  to  heating.  The  term  in  the  square  brackets  of  Eq.  (1)  corresponds  to  the  fraction 
of  the  total  absorbed  photons  that  lead  to  an  excitation  of  the  cooling  transition  in  the  rare- 
earth  ion  (e.g.  2F7/2  -  2F5/2  transition  in  Yb3+),  and  therefore  is  referred  to  as  absorption 
efficiency  t]abs.  The  product  i]extTjabs  can  therefore  be  interpreted  as  the  efficiency  (or 
probability)  of  converting  an  absorbed  photon  to  an  escaped  fluorescence  photon.  The 
background  absorption  ab  most  likely  originates  from  unwanted  contamination  during  growth, 
particularly  from  transition  metals  such  as  copper  and  iron,  similar  to  glass  hosts  [14].  In  the 
cooling  model,  ab  is  also  taken  to  be  temperature  independent,  and  broadband  (i.e. 
independent  of  wavelength)  within  the  spectral  region  of  the  cooling  transition  [14]. 

Due  to  Boltzmann  distribution  of  populations  in  both  ground  and  excited  states,  as  well  as 
the  vibrationally  broadened  levels  in  the  Stark  manifolds,  both  kf  and  a  are  functions  of  lattice 
temperature.  Upon  lowering  of  temperature,  the  cooling  efficiency  in  the  cooling  region  (for 
fixed  k  and  ab )  decreases  due  to  diminishing  resonant  absorption  a  and  red-shift  of  kf  until  it 
becomes  negative,  corresponding  to  overall  heating.  The  temperature  corresponding  to  the 
transition  point  from  cooling  to  heating  is  the  minimum  achievable  temperature  (MAT)  at  that 
wavelength,  i.e.  MAT(l).  For  a  given  material,  a  minimum  temperature  Tmm  of  the  MAT 
spectrum  (i.e.  a  global  MAT)  is  defined  as  the  minimum  temperature  Tmin  at  which  tjc(k,TmirJ  = 
0  or  simply  the  temperature  at  which  dMAT(k)  / dk  =  0  .  Prediction  of  MAT (k)  and  hence 
Tmin  of  a  given  material  requires  knowledge  of  r\c(k,T)  which  in  turn  necessitates  that  its  four 
constituents  a(k,T),  kf  (T),  rjexl  and  ab  be  known.  As  will  be  detailed  in  Section  4,  we 
experimentally  determined  these  four  parameters  for  a  5%  Yb:YLF  crystal  for  the  excitation 
in  the  E||c  orientation.  Figure  1  shows  thus  constructed  temperature -dependent  spectra  of  the 
cooling  efficiency.  A  minimum  achievable  temperature  (MAT)  of  ~115  K  is  predicted  to 
occur  at  -1020  nm,  corresponding  to  the  transition  between  E4  and  E5  levels  of  the  2F7/2  - 
2F5/2  Stark  manifold  of  Yb3+. 

a) 


2r 
r  7Z2 


"1000  1010  1020  1030  1040  1050  1060  1070 

wavelength  (nm) 


Fig.  1.  (a)  Schematic  of  the  Stark  manifold  associated  with  2F7/2  -  2F5/2  transition  in  Yb3+  ion; 
(b)  Calculated  temperature-dependent  spectra  of  the  cooling  efficiency  r\c  of  a  5%  Yb3+:YLF4 
crystal  for  E||c  orientation  with  r\a,  =  0.995  ±  0.001  and  ab  =  (4.0  ±  0.2)  x  lO^1  cnf1,  as 
determined  from  the  fit  of  the  experimental  data  at  300  K  (black  circles).  Transition  from 
cooling  to  heating-only  occurs  between  120  K  (blue  line)  and  1 10  K  (red  line),  bracketing  a 
minimum  achievable  temperature  at  ~1 15  K  near  1020  nm,  corresponding  to  E4-E5  Stark 
manifold  transition  of  Yb31.  Resonant  feature  corresponding  to  E3-E5  transition  is  also  evident 
at  -1010  nm,  but  does  not  lead  to  cooling  efficiency  enhancement  for  given  values  of  r]ext  and 
ab. 
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3.  Two-band  differential  spectral  metrology  (2B-DSM) 

In  order  to  directly  verify  the  cooling  efficiency  spectra  (Fig.  1)  and  in  particular  MAT (/'),  we 
set  out  to  perform  a  pump-probe  experiment  to  measure  “local”  laser-induced  temperature 
change  as  a  function  of  global  (bulk)  temperature  of  the  sample.  For  this,  the  bulk  temperature 
of  the  Yb:YLF  is  varied  by  means  of  a  cryostat  through  an  efficient  thermal  contact  of  the 
crystal  with  the  cryostat’s  coldfinger.  In  this  arrangement,  laser  cooling  occurs  only  within  a 
limited  spatio-temporal  (i.e.  local)  window.  That  is,  the  laser-irradiated  spot  initially  cools 
before  reversing,  as  the  heat  from  fluorescence  absorption  at  the  surface  contacts  diffuses 
back  into  the  cooled  region.  While  the  spatial  variation  of  the  temperature  field  has  been 
addressed  previously  [15,16],  the  temporal  dependence  pointed  out  here  is  due  to  the 
particular  geometry  of  the  experimental  setup,  i.e.  the  presence  of  good  thermal  contact  of  the 
crystal  with  the  coldfinger.  The  duration  of  this  local  cooling  becomes  shorter  for  crystals 
with  relatively  high  thermal  conductivity  such  as  YLF  (as  compared  with  glasses).  This 
restriction  motivated  us  to  develop  a  novel  and  highly  sensitive  fast  fluorescence-based 
temperature  measurement  technique.  This  technique  is  described  in  the  current  section  while 
its  application  to  determination  of  MAT(l)  is  discussed  in  Section  4  of  the  paper. 

Before  we  describe  our  technique,  we  also  should  note  that  a  pump-probe  arrangement  to 
detect  local  temperature  (i.e.  spatio-temporally  resolved  temperature  in  the  region  interrogated 
by  the  probe  beam)  through  an  induced  refractive  index  change  (photo-thermal  deflection 
spectroscopy  PTDS)  was  employed  in  the  past  to  perform  room  temperature  [5,13,17,18]  as 
well  as  cryogenic  [13,17]  cooling  measurements  in  Yb-glasses.  However,  techniques  relying 
on  a  refractive  index  measurement  (e.g.  PTDS)  need  to  separate  the  thermal  contribution  to 
the  index  change  from  the  accompanying  contribution  caused  by  pump-induced  population 
changes  in  both  ground  and  excited  states  [19].  We  note  that  such  separation,  in  principle, 
cannot  be  achieved  by  wavelength  selection  of  the  probe  alone,  since  pump-induced 
absorption  change  leads  to  a  refractive  index  change  at  all  wavelengths  as  per  Kramers- 
Kronig  relations.  This  complication  in  turn  necessitates  a  lateral  displacement  of  the  pump 
and  probe  beams,  since  heat  diffuses  outward  while  the  population  excitation  remains  local. 
Thus  the  need  for  optimal  focusing  geometries  for  both  pump  and  probe  beams  and  their 
optimal  lateral  displacement  (in  general  both  are  temperature  dependent  through  thermal 
diffusivity  of  a  medium,  e.g.  YLF  [20])  complicate  the  measurements  and  can  even  lead  to 
erroneous  results,  especially  at  low  temperatures. 

The  fluorescence  spectrum  of  the  solid  can  also  be  used  to  infer  its  temperature  with  high 
sensitivity  [11,21-23].  In  this  situation,  carrier  density  effects  are  alleviated  by  maintaining 
probe  beam  at  constant  power  and  therefore  no  displacement  between  pump  and  probe  beams 
is  required.  Sensitivity  is  also  increased  using  difference  luminescence  thermometry  (DLT), 
obtained  upon  subtraction  of  the  acquired  spectra  at  two  different  temperatures  [11,22,23]. 
While  highly  sensitive,  the  speed  of  this  technique  is  limited  by  the  spectrometer  (CCD) 
readout  time,  which  is  typically  in  the  tens  of  milliseconds  range  for  the  standard  devices.  In 
addition,  sensitivity  of  the  measurement  is  degraded  from  inherent  CCD  dark  counts  and 
readout  noise. 

To  circumvent  these  limitations,  we  perform  a  dual-band  spectral  differencing  via  a 
balanced  amplified  photodetector  (BAPD,  for  e.g.  a  bi-cell  PD)  directly  in  the  spectral  domain 
[24].  Due  to  the  large  electronic  bandwidth  and  high  common-mode-rejection  of  these 
detectors,  a  fast  and  shot-noise-limited  performance  can  be  obtained  relatively  easily. 
Fluorescence  spectrum  can  be  sampled  for  example  by  means  of  a  monochromator  (MC)  [24], 
color  filters  [25],  or  any  other  spectrally  selective  components.  Color  filters  for  example  are 
more  suitable  for  resolving  a  temperature-induced  spectral  linewidth  change  of  a  fixed 
transition  frequency  (e.g.  Yb3 1  ion),  while  the  MC  offers  a  continuously-tunable  choice  for  the 
temperature-dependent  transition  frequencies,  such  as  bandgap  of  a  semiconductor  (e.g. 
GaAs).  We  have  termed  this  technique  2-band  differential  spectral  metrology  (2B-DSM). 
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In  this  work  we  exploit  the  temperature -dependent  bandgap  shift  of  a  GaAs/GalnP  double 
heterostructure  (GaAs  DHS),  attached  to  the  Yb:YLF  sample  [12]  as  a  means  of  the  2B-DSM 
probe.  To  this  end,  a  portion  of  the  luminescence  spectrum  of  GaAs  S(l)  is  interrogated  at  the 
exit  slit  of  a  monochromator  by  the  B  APD,  such  that  the  resulting  signal  (V)  is  [24] : 


V(A0,ro)XS(A0+AA,ro)-S(A0 
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where  T0  is  the  GaAs  temperature,  k0  is  the  center  wavelength  of  the  spectral  region  within  the 
exit  slit  of  width  2AX.  The  spectral  derivative  approximation  assumes  that  2 A),  is  still  narrower 
than  any  spectral  feature  of  the  luminescence  lineshape  S(A).  For  small  temperature  change 
(AT)  about  the  initial  To,  the  V  signal  (up  to  the  first  order)  is  then  given  by: 


V(AT,A0) 


'  8S 

d 

'  8S 

\ 

y8A 

4»W 

8T 

k) 

AJ+o[(Ar)2]. 


(3) 


By  the  appropriate  choice  of  k0,  temperature-independent  offset  [first  term  on  the  right 
hand  side  of  Eq.  (3)]  can  be  nullified,  such  that  the  signal  V  (in  both  magnitude  and  sign)  is 
background-free  and  linearly  proportional  to  the  temperature  change  AT.  We  note  that  even  if 
the  derivative  approximation  in  above  expressions  cannot  be  made,  spectral  shift  of  the 
luminescence  as  measured  by  the  differential  V  signal  can  still  be  calibrated  to  be  proportional 
to  the  temperature  change. 


Fig.  2.  (a)  Schematic  of  the  pump-probe  setup,  BAPD  -  balanced  amplified  photodiode,  signal 
(V)  is  proportional  to  the  change  of  temperature  (AT);  (b)  Normalized  calculated  and  box-car 
averaged  thenno- spectral  derivative  (orange  line)  of  a  GaAs  DHS  evaluated  at  Ao  and  To  = 
300K  is  plotted  versus  Ao  together  with  the  normalized  magnitude  of  a  measured  differential 
modulated  signal  (V)  (black  circles);  (c)  Modulated  pump  (green)  and  temperature  (red),  as 
obtained  from  the  probe  signal  when  GaAs  semiconductor  is  excited  by  50  mW  of  532nm 
pump  and  monitored  at  Ao  ~870  nm;  (d)  Magnification  of  the  panel  (c)  around  t  =  0  point; 
resolution  of  ~250  pK  is  demonstrated  after  10,000  waveform  averages  on  an  oscilloscope. 


First  we  perform  stand-alone  experiment  on  GaAs/InGaP  double  heterostructure  to 
demonstrate  the  sensitivity  of  the  2B-DSM  technique.  Here,  a  500  pm  thick  GaAs  DHS  is 
excited  from  opposite  sides  by  a  pump  laser  (50  mW  at  532  nm,  modulated  at  10  Hz)  and  a 
650  nm  probe  laser,  focused  to  their  respective  spot  sizes  of  ~200pm  and  150pm  [Fig.  2(a)]. 
Probe-induced  luminescence  is  collected  by  a  600  pm  diameter  optical  fiber  and  routed  into 
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the  entrance  slit  of  a  scanning  MC.  Closely  spaced  and  properly  oriented  cores  of  two  optical 
fibers  at  the  exit  slit  of  the  MC  pick  up  the  respective  components  of  the  spectrally-dispersed 
luminescence  signal  and  couple  them  into  the  input  ports  of  the  BAPD  (New  Focus  Model 
2307).  Modulated  pump  induces  a  time -varying  differential  signal,  which  is  pre-amplified  and 
then  averaged  on  an  oscilloscope  [Fig.  2(c)].  Complimentary  to  the  differential  signals, 
spectrometer  (Ocean  Optics  FIR4000  with  10pm  slits)  is  used  to  measure  luminescence 
spectra  of  GaAs  DHS  with  and  without  modulated  pump.  These  additional  measurements 
allowed  us  to  calibrate  pump-induced  heating  of  the  sample  (through  a  Varshni  relation  [21]) 
as  well  as  allowed  for  a  comparison  of  the  differential  signals  with  that  of  the  predicted 
behavior  in  Eq.  (3).  Figure  2(b)  shows  room-temperature  calculated  spectrum 

of  d2S  /  dA8T |  ,  performed  on  a  box-car  averaged  spectra  S(X),  with  a  box  size  corresponding 

to  Al  of  5  nm  of  the  MC.  Maximum  signal  deviation  occurs  for  hv„  ~Eg,  the  bandgap  energy 
of  the  GaAs  DFIS  (10  =  c/vo  -870  nm).  Magnitude  of  the  time -dependent  component  of  the 
differential  signals  is  plotted  on  the  same  axes  for  various  grating  positions  (10).  A  very  good 
agreement  is  observed,  verifying  Eq.  (3)  and  in  particular  the  linear  dependence  of  the 
differential  signal  (V)  on  the  small  AT,  induced  by  the  pump  beam.  An  uncertainty  in  abscissa 
coordinate  of  the  2B-DSM  signal  is  mainly  due  to  the  backlash  of  the  screw-controlled  grating 
position  of  our  compact  mini-MC  (Optometries).  Small  systematic  differences  between 
calculated  and  measured  values  of  the  derivative  spectra  are  attributed  to  spectral  response 
differences  of  the  spectrometer  CCD  and  each  of  the  balanced  detectors,  together  with 
possible  higher  order  contributions  (neglected  terms  in  Eq.  (3). 

To  identify  the  temperature  resolution  of  the  technique,  we  average  10,000  time -resolved 
oscilloscope  traces  of  the  probe  signal  monitored  atl0=  870  nm,  and  converted  to  temperature 
by  means  of  a  Varshni  formula  [21].  For  a  50  mW  pump  laser,  heating  of  the  heat-sunk  GaAs 
DHS  by  0.6  degrees  was  deduced  [Fig.  2(c)].  A  close-up  view  at  the  early-time  probe  signal 
reveals  temperature  resolution  of  -250  pK  over  on  a  sub-millisecond  time  scale  [Fig.  2(d)]. 
This  performance  allows  for  a  sensitive  determination  of  the  MAT  {'/.)  spectrum. 

4.  Spectroscopy  of  the  minimum  achievable  temperature  (MAT) 

The  spectroscopy  of  the  cooling  efficiency  is  performed  on  a  high-purity,  Czochralski-grown 
5%-doped  Yb:YLF  crystal  of  dimensions  3  x  3  x  9  mm3  [11,12,26],  excited  in  the  E\\c 
orientation  along  the  long  axis.  As  was  mentioned  above,  for  modeling  of  the  cooling 
efficiency  (Fig.  1),  we  experimentally  obtained  quantities  a.(l,  T),  If  (T),  t]exl  and  ab.  To  this 
end,  temperature -dependent  polarized  and  unpolarized  instrument-response-corrected 
fluorescence  spectra  of  the  Yb:YLF  sample  were  obtained  by  placing  the  crystal  inside  a 
cryostat.  Mean  fluorescence  wavelength  l/T)  was  calculated  by  taking  a  first  moment  of  the 
unpolarized  fluorescence  spectra,  while  resonant  absorption  spectra  a(X,T)  were  obtained  from 
the  reciprocity  analysis  of  the  polarized  emission. 

To  obtain  external  quantum  efficiency  and  background  absorption  values,  an  additional 
experiment  was  performed  at  room  temperature,  where  sample  was  excited  by  a  tunable 
(950nm-1080nm,  1.3-1. 8  W)  CW  Ti:Sapphire  laser  and  crystal  temperature  was  monitored 
using  a  bolometric  thermal  camera.  The  obtained  temperature  change,  normalized  by  the 
absorbed  power,  is  linearly  proportional  to  the  cooling  efficiency.  We  fit  the  data  (Fig.  1, 
black  circles)  to  the  cooling  efficiency  relation  [Eq.  (1)]  to  extract  tjext  =  0.995  ±  0.001  and  ab 
=  (4.0  ±0.2)  x  10-4  cm-1. 

Next,  we  discuss  the  application  of  a  2B-DSM  technique  to  test  modeling  predictions  in 
Fig.  1.  The  tunable  pump  (Ti:Sapphire  laser)  excites  Yb:YLF  sample  held  at  variable  starting 
temperatures  (To)  by  means  of  a  contact  with  a  coldfinger  of  a  liquid  nitrogen  filled  cryostat 
[Fig.  3(a)],  which  uses  a  thermal  impedance  displacer  together  with  a  built-in  heating  element 
to  allow  for  a  potential  set  temperature  in  the  range  of  80  -  350  K.  Luminescence  of 
GaAs/GalnP  DHS,  attached  to  the  side  of  the  Yb:YLF,  is  excited  by  a  CW  probe  beam  (650 
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inn,  20  mW)  and  collimated  by  a  lens  (N.A.  =  0.42)  after  the  cryostat  viewport.  Another  lens 
(N.A.  =  0.22)  is  used  to  refocus  the  luminescence  signal  into  a  NA-matched  optical  fiber  (600 
pm  core),  which  in  turn  is  coupled  into  an  entrance  port  of  a  mini-monochromator  for 
differential  detection  (as  detailed  in  the  previous  section). 

Pump-induced  temperature  change  AT,  corresponding  to  laser  cooling,  is  expected  to  be 
small  and  fast,  following  initial  pump  turn  on.  This  is  because  at  times  larger  than  the  thermal 
transport  time  across  the  crystal,  its  temperature  is  dictated  not  by  the  pump  laser  but  by  the 
heat  generated  by  fluorescence  absorption  at  the  coldfinger,  diffusing  from  the  coldfinger 
contacts  back  into  the  excitation  volume.  Fast  measurement  of  the  DHS  temperature  is  thus 
capable  of  separating  the  initial  local  temperature  T0  -  AT  (due  to  laser  cooling  process)  from 
the  background  heat  at  the  later  times.  An  additional  thermal  capacitance  has  been  added  at 
the  contact  points  of  the  crystal  with  the  coldfinger  surfaces  (150  pm  thick  BK7  plates)  in 
order  to  increase  the  transport  time  of  the  background  heat  and  thus  increase  the  time  window 
over  which  the  local  temperature  change  can  be  observed.  Finally,  pump  is  traversing  the 
Yb:YLF  parallel  and  very  near  the  surface  with  the  attached  DFIS  [Fig.  3(a)],  thus  minimizing 
thermal  transport  time  of  the  laser  cooling  signal,  as  detected  by  the  bandgap  shift  of  the  DFIS. 

Figure  3(b)  shows  time  traces  of  the  2B-DSM  signals  (spectral  derivative,  see  Section  3) 
from  Yb:YLF/GaAs  at  room  temperature.  When  pump  is  tuned  below  (980  nm)  and  above 
(1020  nm)  the  mean  luminescence  wavelength  (-1000  mn),  clear  phase-shift  of  the  time  traces 
is  evident.  For  a  given  excitation  wavelength,  the  transition  temperature  between  local  cooling 
and  heating  determines  the  minimum  achievable  temperature  (MAT)  for  that  excitation 
condition.  For  each  excitation  wavelength,  we  vary  the  bulk  sample  temperature  ( T0 )  via  the 
cryostat,  until  transition  from  cooling  to  heating  is  observed,  corresponding  to  a  sign  change 
of  the  measured  AT  [Fig.  3(b)].  These  transition  temperatures  are  then  plotted  in  Fig.  3(c), 
superimposed  on  top  of  a  contour  plot  of  the  cooling  efficiency  spectra,  as  calculated  from  Eq. 
(1)  and  also  represented  in  Fig.  1.  The  red  color  corresponds  to  heating  regions,  while  the  blue 
area  is  consistent  with  cooling.  The  line  separating  these  regions  defines  the  MAT(1). 


time  (ms)  wavelength  (nm) 


Fig.  3.  (a)  Schematic  of  the  experimental  arrangement:  Yb:YLF  is  clamped  by  the  coldfinger 
arrangement  that  is  held  at  To;  local  temperature  change  due  to  the  pump  beam  is  detected  via 
luminescence  from  a  GaAs/InGaP  double  heterostructure,  excited  in  turn  by  the  probe  laser. 
Thermal  buffers  serve  to  maximize  local  signal,  while  maintaining  Yb:YLF  temperature  near 
the  To  setpoint;  (b)  Normalized  and  vertically-shifted  time  traces  of  the  spectral  derivative 
signals,  showing  a  clear  phase-reversal  between  the  heating  (980  nm)  and  cooling  (1020  nm) 
excitations  at  room  temperature;  (c)  comparison  of  the  contour  plot  of  the  calculated  cooling 
efficiency  with  the  measurement  (circles)  of  the  wavelength-dependent  minimum  achievable 
temperature  MAT(^),  separating  cooling  (blue)  and  heating  (red)  regions,  cooling  to  1 10  ±  5  K 
at  ~1020  nm  is  demonstrated;  inset  shows  energy  diagram  (not  to  scale)  of  the  Yb3+  Stark  level 
manifold. 


We  point  out  that  MAT  condition  predicted  from  the  model  is  in  a  very  good  quantitative 
agreement  with  the  measurements  [Fig.  3(c)],  verifying  both  the  narrowing  of  cooling  window 
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at  lower  temperatures  and  its  convergence  to  the  wavelength  of  1020  nm,  corresponding  to  the 
E4-E5  Stark  manifold  transition  [11].  In  particular,  a  measured  MAT(1020  nm)  =  1 10  ±  5  K  is 
in  excellent  agreement  with  the  model,  within  the  current  experimental  uncertainty,  as 
determined  by  the  experimenters’  choice  of  the  temperature  step  in  T0.  This  temperature  is 
below  the  NIST -defined  cryogenic  onset  of  123  K  and  corresponds  to  the  lowest  temperature 
obtained  by  means  of  optical  refrigeration.  Other  researchers  have  also  reported  or  inferred 
observation  of  local  cooling  below  110  K  in  ytterbium-doped  glasses  [13,17,18]  using  photo- 
thermal  deflection  spectroscopy.  As  discussed  above,  however,  photo-thermal  deflection 
measurements  are  difficult  to  interpret  unambiguously,  especially  near  the  MAT  where  the 
signal  to  noise  is  small.  The  present  results  are  not  plagued  by  these  ambiguities.  In  particular 
we  have  demonstrated  quantitative  agreement  between  the  empirical  modeling,  local 
measurements  and  the  earlier  reported  bulk  cooling  of  Yb:YLF  to  155  K  at  1023  nm  [11]. 
Demonstrated  success  of  the  cooling  efficiency  model  also  allows  us  to  more  confidently 
estimate  MAT  of  a  typical  Yb:ZBLAN  (fluorozirconate  glass  host)  to  be  — 180  K.  This  value 
is  warmer  than  MAT  of  Yb:YLF  due  to  large  inhomogeneous  broadening  of  the  Stark- 
manifold  levels  in  glass  host,  together  with  smaller  allowable  Yb3+  doping  concentration  (than 
in  YLF)  and  comparable  background  absorption  values  [27].  Furthermore,  previous  bulk 
cooling  measurements  in  Yb:ZBLAN  to  208  K  at  1026  nm  [27]  are  also  consistent  with  the 
estimates  from  the  current  cooling  efficiency  model. 

Finally,  we  note  that  this  verification  of  the  model  also  supports  the  built-in  assumptions, 
namely  wavelength-  and  temperature-independence  of  the  background  absorption  coefficient. 
This,  in  turn,  suggests  transition  metal  impurities  as  the  likely  source  of  the  parasitic 
background  absorption  in  YLF,  in  similarity  to  the  glass  hosts  [14].  Identifying  these 
impurities  directly  is  a  formidable  task,  as  their  estimated  concentrations  are  well  below  part 
per  million  levels.  The  obvious  incentive  however  is  that  temperatures  below  100  K  and  even 
approaching  77  K  are  possible  upon  an  order-of-magnitude  improvement  of  the  present 
crystal  purity  [28]. 

5.  Conclusion 

Measurement  of  a  spectrum  of  the  minimum  achievable  temperatures  in  a  5%-doped  Yb:YLF 
crystal  is  reported  using  a  novel  pump-probe  technique.  The  obtained  spectrum  and  in 
particular  its  minimum  of  110  ±  5  K  at  E4-E5  Stark  manifold  resonance  of  Yb3+,  are  in 
excellent  agreement  with  the  model  predictions.  This  temperature  corresponds  to  the  coldest 
(local)  temperature  obtained  by  optical  refrigeration  to  date  and  is  below  the  NIST-defined 
onset  of  cryogenics  at  123  K. 

These  results  were  possible  via  a  non-contact  and  highly  sensitive  2-band  differential 
spectral  metrology  (2B-DSM,  patent  pending)  pump-probe  technique.  In  particular  we  have 
demonstrated  temperature  resolution  of  ~250  pK  on  a  sub-millisecond  time  scale.  The  speed 
and  sensitivity  of  this  technique  have  allowed  us  to  resolve  the  initial  (local)  temperature 
dynamics  in  the  Yb:YLF  crystal.  This  technique  is  general  and  we  envision  applications  in 
derivative  spectroscopic  methods  as  well  as  flow  cytometry  to  name  few  examples. 

While  this  work  verifies  the  predicted  potential  of  the  current  Yb:YLF  sample  to  cool  to 
— 110  K,  bulk  cooling  from  ambient  to  this  temperature  using  a  high  power  laser  at  -1020  nm 
remains  to  be  demonstrated. 
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